Second harmonic generation in multilayer graphene induced by direct electric current 
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Optical second harmonic generation (SHG) is studied from multilayer graphene films in the pres- 
ence of DC electric current flowing in the sample plane. Graphene layers are manufactured by 
chemical vapour deposition (CVD) technique and deposited on an oxidised Si(001) substrate. SHG 
intensity from graphene layer is found to be negligible in the absence of the DC current, while it 
increases dramatically with the application of the electric current. The current-induced change of 
the SHG intensity rises linearly with the current amplitude and changes its sign under the reversal 
of the current direction to the opposite. The observed effect is explained in terms of the interference 
of second harmonic radiation reflected from the Si surface and that induced by the DC current in 
multilayer graphene. 
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Since its first experimental realisation in 2004 graphene continues to attract enhanced interest as a prospective 
material for both fundamental and applied science. Fascinating electronic properties which include electric field- 
effect pi, "chiral" quantum Hall effects @, Q) prospects for spintronics [J] and valeytronics Q immediately pushed 
graphene research to the cutting edge of modern nanomaterial science and technology. Among the numerous problems 
currently being studied for graphene is the possible connection between the electron transport and the nonlinear-optical 
response. The importance of this task is dictated not only by needs of the applied research as allows distant probing 
of the electron flow in graphene devices but, perhaps, more importantly as a route to gain new comprehensive insight 
into its fundamental electronic properties. 

Second harmonic generation (SHG) is among the most ubiquitous methods used for probing surfaces and interfaces of 

centrosymmetric materials. [6| High sensitivity to the surface and thin film properties arises from SHG being prohibited 

in the electric dipole approximation in the volume of a centrosymmetric medium. As a result it is generated basically 

£h ■ at surfaces and interfaces where the central symmetry is broken. Moreover one can break the inversion symmetry by an 

external influence such as electric and magnetic field causing so-called field induced second harmonic generation. [7H9| 

i^i . It has been demonstrated recently both theoretically and experimentally[10|, LU| that DC electric current flowing in 

the plane of a centrosymetric semiconductor can break the symmetry of the electron density distribution, resulting in 

t~H , current-induced SHG (CSHG) which can overwhelm conventional electric-field-induced mechanism if the conductivity 

k* - ■ of the probed material is sufficiently high. Moreover, theoretical predictions [ll| made almost a decade before the 

advent of graphene demonstrate the possibility of the SHG enhancement by 1~2 orders of magnitude in case of 

ji- ■ ballistic electron transport and in case of two-dimensional nature of the investigated electron system. In this paper 

_o- ' we report the first investigation of current-induced second harmonic generation in multilayer graphene under ambient 
• , conditions. 

^-J ' Our specific experimental conditions are the following. Either p, s -polarised output of a femtosecond ThSapphirc 

laser system operating in the wavelength range of 730^830 nm, pulse duration 100 fs, at 80 MHz repetition rate 

,_^ ■ is focused into a 50 fim size spot on the sample at 45° angle of incidence with the peak intensity of approximately 
0, 8GW/cm 2 . SHG radiation reflected in the direction of specular reflection was spectrally selected by BG39 color 
Schott filters and detected by a PMT. Sample holder could be rotated around the axis orthogonal to its surface thus 
allowing the azimuthal SHG studies. Application of the DC electric current was performed by mechanical adjusting of 
two spring-assisted Pt needles to the graphene film distanced by 2-^3 mm. To ensure the the second order nonlinearity 
is the source of the signal at the SH frequency the conventional quadratic dependendence of the measured signal on 
the fundamental radiation intensity was checked. 

Multilayer graphene film of 4^5 monolayers thickness was manufactured by chemical vapour deposition technique 
on a thin polycrystalline Ni foil and transferred to SiC>2/Si(001) substrate with a 300 nm oxide layer fl2Hl4j |. The 
size of the graphene-coated area was 7x7 mm 2 approximately. The size of single crystal domains in the film was 
shown to be in the range of 3-5 fim as evaluated from the SEM data and that is typical for this method of graphene 
composition |13j. Such graphene films possess high electrical conductivity and is capable of preserving linear electron 
dispersion law [15l4l7j ] due to stacking disorder. 
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FIG. 1: (Color online) Azimuthal dependencies of the SHG intensity for bare Si(001) (green (light gray) line) and for graphene- 
coated Si(001) (black line); Fourfold sine fit (blue (dark gray) dashed line). DC electric current J—0 
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FIG. 2: Current dependence of the SHG contrast obtained for pp combination of polarisations of the fundamental and SHG 
waves for A =800 nm. Inset: Time dependence of the SHG intensity, CSHG contrast for inset « 0.70 



As SHG is sensitive to the symmetry of the probed medium it is crucial to define the symmetry of the observed effect. 
Figure [T] shows a typical SHG azimuthal pattern measured in order to define the graphene film crystal anisotropy. 
Fourfold symmetry of the SHG intensity pattern is observed that is governed by the SHG anisotropy from Si(001), 
while graphene film with the structural symmetry 6mm should reveal three- or sixfold symmetric SHG pattern [18|. 
Thus we have to consider that the main nonlinear-optical source of the sample is the Si(001) surface. 

A steep rise in the SHG intensity is observed as the DC current is applied (Fig. [51 inset). To characterise the 
relative magnitude of the effect we used current-induced SHG contrast, p = 2(/ 2w (J) — 7 2 "(0))// 2 "(0), where 7 2 "(0) 
and I 2ul {J) are the SHG intensity measured in the absence and in the presence of the current, J. Figure [5] shows that 
p( J) dependence is linear within the experimental accuracy. This effect is consistent with the model SHG description 
based on the interference of the two SHG fields, I 2uJ oc (E 2uj + E 2u (J)) 2 « const + 2E^. aphene (f)E^Cos(A<f>) oc J, 
where Eq t aphene and Egf are electric fields of second harmonic waves from graphene ans silicon, respectively, and 
Acf> is the phase shift between these fields, J is the current value. Here we assume that current-induced SHG is much 
lower as compared with that generated by Si(001) surface. 

The symmetry of the CSHG can be estimated for a material with any given crystal symmetry without even 
specifying the microscopic mechanism that underlines its inception. CSHG electric-dipole polarisation is given by 
Pf u = Xijki EfEfr J;, where E is the fundamental electric field and Xlju is an effective four-rank susceptibility 
tensor which governs the SHG effect. It can be easily shown (see appendix) that the dependencies of the CSHG 
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FIG. 3: CSHG contrast dependencies on the current value for different combinations of polarisations of the pump and SH waves 
(A=800 nm) for (a) 0=0° and (b) </>=90°; Lines are linear fit of the experimental data 

intensity on the angle </> between the direction of the current flow and the plane of incidence are I PPlS p °c acos0| 2 , 
and I ps ,ss oc |6sin0| 2 for pp,sp and ps,ss combinations of the input and SHG polarizations, respectively. Here a,b 
are the combinations of the j^ 2 ^ components. Thus the CSHG contrast should depend on the polarization of the 
pump and SHG beams. 

Figure [3] shows the dependencies of the CSHG contrast on the current for different <j> values. It can be seen that 
in case of the current flow being parallel to the plane of incidence (0=0°) the CSHG contrast for p-polarizcd SHG 
is much higher as compared with s-polarized one, while for 0=90° the situation is inverted. This stays in agreement 
with the symmetry description discussed briefly above. 

Finally we performed the SHG spectroscopy measurements. Figure U shows the SHG intensity spectra obtained 
in the absence of the DC current as well as its current-induced modification. In the former case the SHG spectrum 
is consistent with that of a Si(001) surface. It demonstrates a maximum centered at approximately 3,34 eV (A k, 
740 nm) which corresponds to the two-photon direct interband transition in silicon J19j. The SHG spectrum shifts 
significantly under the application of the DC current with the sign of this shift depending on the direction of the 
current . 
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FIG. 4: (Color online) (a) SHG intensity spectrum: Black dots - SHG intensity spectrum without current, red and blue dots 
- SHG intensity spectrum with current in different directions (b) CSHG electric field (open squares) and phase shift (closed 
circles) spectra; Inset: Transmission spectra of the mutilayer graphene on glass in the near-UV 



These spectral dependencies can be understood considering a phase shift between the interfering SH signals from 
the silicon surface and CSHG from graphene which varies with the wavelength. To prove this explanation right 
we have calculated the CSHG electric field and phase shift spectra from fig. [4^ assuming that the phase of the 
current-induced SHG wave changes by tt with the change of the direction of the current to the opposite. The phase 
shift (Fig. lib, circles) resembles conventional spectral shape in the vicinity of the resonance in silicon similar to 



what we have expected. The more remarkable however is the pure CSHG spectrum (Fig. [4}j, open squares) which 
shows an enhancement at the short wavelength edge correlating with the rise of the absorption at the corresponding 
two-photon energy (Fig. HJd, inset). This is generally believed to come from high energy excitonic resonances at 4,5 
eV [20, [21| and/or interband transitions at 5,1 eV[22|. Thus we believe the same mechanisms to play their role in 
the CSHG enhancement in graphene. Still with no theory proposed we cannot specify the microscopic mechanism of 
such enhancement and this subject seeks further investigation. 

It is also necessary to mention that the observed effect may be the result of electric-field induced SHG (EFISH) driven 
by planar field associated with the current flow as a remarkably high \ has been recently reported for graphenc[23]. 
However we believe that CSHG should overwhelm the EFISH in graphene because of high conductivity(resulting in 
insufficiently high electric field). 

In conclusion we have measured second harmonic generation from multilayer CVD graphene subjected to electron 
flowing within the sample plane. As the film itself was found to produce no significant SH response the situation 
changed drastically with in plane current application. The dependence of the CSHG effect magnitude on current 
density is found to be linear in a good agreement with the theory and the effect was observed to depend on current 
directions as was expected from the symmetry analysis. These two results combined provide a clear evidence of it 
being possible to use the discussed technique for distant probing of current density distribution in graphene devices. 
Finally a significant interference-mediated SHG spectral shift associated with current flowing in the sample and the 
CSHG enhancement at the short wavelength spectral edge possibly connected to resonances in graphene in UV are 
observed. 
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Appendix: CSHG symmetry analysis 

The symmetry of the CSHG effect it determined by the symmetry of the effective forth rank susceptibility tensor 
X^ 2,1 ) that governs the current-induced component of the SHG field in the electric-dipole approximation, E 2u (J) oc 
P 2uj (J) — x ' ] E U E" J . While an ideal graphene layer possesses 3m symmetry, it stems from our SHG azimuthal 
studies that no anisotropic SHG comes from the graphene layer in the absence of DC current. This may be due to 
the fact that the SHG measured in the experiment is the result of averaging over the laser spot (of about 50 /Ltm 
in diameter) which contains hundreds of chaotically oriented crystallites each of them being several microns in size. 
Thus we can consider the symmetry of the graphene layer as that of an isotropic surface. 

As the second step we will use the same symmetry of the CSHG and EFISH (electric field induced SHG) [HJ and 
this find nonzero x tensor elements for the chosen experimental geometry. An important fact here is that these 
CSHG components are odd in J, i.e. change their sign under the reversal of the current direction to the opposite. Two 
coordinate system that are considered below are shown in Fig. [5j the laboratory one (XYZ) and that one connected 
with the sample, (X'Y'Z'). The axes Z and Z' are oriented along the normal to the surface while (XOY) and (X'OY') 
are parallel to the surface of the sample. Azimuthal rotation of the sample results in the modification of the angle <f> 
between the X and X' axes. The plane of incidence corresponds to (XOZ), p and s polarizations of the light waves 
are also shown in the Figure. 
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FIG. 5: Surface and laboratory coordinate systems 

In order to reveal the azimuthal CSHG dependence we first consider nonzero X*- 2 ' 1 -* components for the DC current 
J ||X',X (0=0), which are Xz>z'x>x>,Xz>x'z'x',Xx'z>z>x',Xx'x'x'x' for the p-in, p-out polarizations. Let us consider the 



azimuthal rotation of the sample in that case. The SHG field in laboratory coordinate system generated by each 

of the x^ 2 ' 1 ' nonzero components is determined by E 2 " = aw djj' dkk' aw ■ Xi'j'k'i'^f'^k'^l'^ where a is a standard 
rotation matrix and the summation over the iterative indices is assumed. This gives us the following expression for 
the rotational anisotropy of the CSHG field Ep SHG ex cos<p. Analogously one can easily obtain E GSHG ex sine/). 

Thus the azimuthal dependence for the CSHG contrast can be estimated. Taking into account that E 2 " = E (J) + 
E 2ul (Q) one can easily show that the CSHG contrast for the p-in, p-out combination of polarizations is determined by 
the ratio 

P PP ,s P = 2(/ 2 "(J) - I 2 "(0))// 2 "(0) ex (A.l) 

E 2ul (J)/E 2u (0)(xcos(t), (A.2) 

P ps ,s8 <x sin<f>, (A. 3) 

In other words, maximal CSHG contrast in pp and sp combinations of polarizations should be observed as the DC 
current is flowing in the plane of incidence, J ||X and should vanish in case J ||Y. Analogously, for the case of 
s-polarized SHG the CSHG contrast should be maximal for J ||Y and vanish if J ||X. 
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